Tunable anisotropic magnetism in trapped two-component Bose gases 
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We theoretically address magnetic ordering at zero and finite temperature in both homogeneous 
and trapped Bose-Bose mixtures in optical lattices. By using Bosonic Dynamical Mean-Field Theory, 
we obtain the phase diagram of the homogeneous two-component Bose-Hubbard model in a three- 
dimensional (3D) cubic lattice, which features competing magnetic order of XY-ferromagnetic and 
anti-ferromagnetic type in addition to the Mott and superfluid states. We show that these magnetic 
phases persist also in the presence of a harmonic trap. 
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I. INTRODUCTION 

Quantum magnetism is one of the most intriguing ar- 
eas in condensed-matter physics. Even though the at- 
tempts to understand it root back to the early days of 
quantum theory, still there are many open questions even 
at the level of the minimalistic Hubbard model [l] . Many 
theoretical and experimental efforts have been devoted 
to revealing the mechanisms behind magnetic ordering 
of many-body systems [2j . Due to the high level of com- 
plexity in solid-state systems, a quantitative comparison 
between theory and experiment seems a very challenging 
task, if not impossible at all. Therefore it is highly desir- 
able to work with systems which are able to simulate the 
original solid-state many-body systems, but in a much 
more controllable way. 

Over the past decade cold-atomic quantum gases in 
optical lattices have provided an excellent laboratory for 
investigating many-body quantum systems with an un- 
precedented level of precision and control. Nowadays 
one can routinely create optical lattices with different ge- 
ometries that mimic solid-state crystalline structures and 
then load different quantum gases into them, playing the 
role of electrons. Based on this experimental progress, 
fundamental many-body phenomena such as the Mott 
insulator quantum phase transition of interacting bosons 
have been realized [3j. For bosonic gases in optical lat- 
tices correlated atom tunneling Q and superexchange 
due to second-order atom tunneling [f| have been ob- 
served, which are the basic mechanisms leading to quan- 
tum magnetism. At the current stage, detecting anti- 
ferromagnetic long-range ordering of spinful fermions or 
bosons is arguably the most challenging goal and there 
are large experimental and theoretical efforts directed to- 
wards reaching this phase. This may eventually give in- 
sights into the mechanism of high-T c superconductivity 
as well The recently developed optical quantum gas 
microscope @, [1] offers the possibility to detect quantum 
magnetism on the single-atom level, and has been uti- 
lized to observe the phase transition of a one-dimensional 
chain of interacting Ising spins by using a Mott insulator 
of spinless bosons in a tilted optical lattice Q . 

In parallel to these experimental developments, there 



are several theoretical proposals how to gain fundamen- 
tal insight into quantum magnetism via ultracold gases 
One of the aims has been to realize a spin Hamil- 
tonian by interacting ultracold bosons [Hi[ll|. For this 
two-component Bose gas in an optical lattice, which is ef- 
fectively described by the Bose-Hubbard model, the com- 
plete phase diagram was mapped out by a variational 
approach pjj. 

A related important experimental achievement is given 
by mixtures of different quantum gases, whether with 
different statistics (Bose-Fermi mixture) or different iso- 
topes of the same type (Bose-Bose or Fermi-Fermi mix- 
ture). Recently a Bose-Bose mixture of 87 Rb and 41 K 
in a 3D optical lattice has been realized [20(. Here it 
is possible to tune independently the inter-species and 
intra-species interactions. In a follow-up experiment the 
same Bose-Bose mixture has been used to investigate 
the entropy exchange between two spin-dependent traps 
though without lattice |2lj. Another Bose-Bose mixture, 
made of two hyperfine states of 87 Rb, has served as a spin 
gradient thermometer which allows to measure the tem- 
perature of ultra-cold atoms in optical lattices [12, HH . 
In a further experiment the effect of the Bose-Bose inter- 
species interaction on the bosonic superfluidity of one of 
the components in a 3D optical lattice has been explored 

Motivated by these experiments, here we theoretically 
investigate a two-component Bose gas in 2D and 3D 
optical lattices. This system can be effectively mod- 
eled by the Bose-Hubbard model. We will specifically 
consider the case of filling n = 1 and n — 2 per 
site. We investigate the homogeneous (untrapped) sys- 
tem by means of Bosonic Dynamical Mean-Field Theory 
(BDMFT) [H-SH and the harmonically trapped case by 
its real-space generalization (RBDMFT), which extends 
the original BDMFT formalism to the study of inho- 
mogeneous systems. DMFT has been established as a 
powerful tool for strongly-correlated fermionic systems 
[3lj . while BDMFT has recently been shown to provide 
a qualitatively and in 3D even quantitatively accurate 
picture of the Bose-Hubbard model [29[ . For the homo- 
geneous system, we extend our previous BDMFT analy- 
sis to higher filling n = 2 and a different set of interac- 
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tion parameters motivated by experiments. We map out 
the phase diagram and obtain diverse phases such as su- 
perfluid, unordered Mott state and XY-ferromagnetic or- 
der. In addition we turn to the inhomogeneous (trapped) 
Bose-Hubbard model which is more closely related to 
the experimental situation. We include the effect of 
the external confining potential by RBDMFT, which as- 
sumes site-dependent self-energies. In parallel we per- 
form also a complementary calculation based on a Local 
Density Approximation (LDA) combined with BDMFT 
which is computationally more affordable. Comparing 
results of both methods, we examine the magnetic prop- 
erties of the system for a wide range of parameters. To 
our best knowledge this is the first systematic and non- 
perturbative study of the magnetic properties of a two- 
component inhomogeneous Bose-Hubbard model. It will 
bring more insight into ongoing experiments on Bose- 
Bose mixtures in optical lattices. 

The paper is organized as follows: in section II we give 
a detailed description of the RBDMFT approach. Section 
III covers our results for the homogeneous Bose-Hubbard 
model. In Section IV we consider the trapped system 
and present results by RBDMFT and LDA+BDMFT. 
We summarize with a discussion in Section V. 



II. MODEL AND METHOD 

We consider a two-component bosonic mixture in an 
optical lattice with either 2D square or 3D cubic geome- 
try. Experimentally the Bose-Bose mixture could consist 
of two different species, e.g. 87 Rb and 41 K as in Ref. 
[20l | or two different hyperfine states of a single species, 
e.g. 87 Rb as in Refs. [1, III III- In addition we include 
an external harmonic trapping potential which gives rise 
to inhomogeneity. This system can be described by a 
two-component inhomogeneous Bose-Hubbard model 

U = - ^] tu{b\ v bj V + h.c) + ^ ^ U\yh it \{h i>v - S X u) 

<i,3> 

v—b^d 

+ ( V i-Vu)fliu (1) 

In this Hamiltonian represent the nearest neighbor 
sites i,j and we denote the two bosonic species as b, d 
which are labeled by the index \{v) — b, d. The bosonic 
creation (annihilation) operator for species v at site i 
is b\ v {bi V ) and the local density is fii M — b\ u bi V . Due 
to possibly different masses or a spin-dependent optical 
lattice, these two species in general hop with non-equal 
amplitudes tb and td- U\ v denotes the inter- and intra- 
species interactions, which can be tuned via a Feshbach 
resonance or by spin-dependent lattices. /j„ denotes the 
global chemical potential for the two bosonic species and 
Vi the harmonic trap. 

For the fermionic Hubbard model, DMFT has been de- 
veloped and implemented successfully (30L l3lj as a non- 



perturbative formalism to study strongly-correlated elec- 
tronic systems. This includes real-space generalizations 
of DMFT to address inhomogeneous fermionic systems 
[32l [331 ] - In spite of this success, a bosonic version of 
DMFT (BDMFT) for the Bose-Hubbard model has been 
formulated [1^] and implemented (2(| [27|, only very 
recently. Inspired by the case of fermions, here we ex- 
tend BDMFT 26] to a real-space BDMFT (RBDMFT) 
formalism for inhomogeneous systems to include also the 
harmonic trap which is crucial in the experiments. RB- 
DMFT is capable of providing an accurate and non- 
perturbative description of the ground-state of the in- 
homogeneous Bose-Hubbard model ([T]). Like fermionic 
DMFT, also BDMFT assumes that the self-energy of the 
system is local. However, in an inhomogeneous system it 
depends on the lattice site, i.e. = where 

Si j is a Kronecker delta. 

In a more formal language, first we map the Hamilto- 
nian ((TJ) onto a set of individual single-site problems each 
of which is described by a local effective action 

Sjg =fdrdr' £ bi%)t^T -r'rMV) 

J ° \M={b,d} 

+ r^{E^^4 i;) (r)(nW(r)-«5 A; ,) 

E^^wv^w+ft^w^crr)} (2) 
<*j>.f 

where i is the index of the lattice site. In this equation 
t is imaginary time and the function G^\ v (j — t') is 
a local non-interacting propagator interpreted as a local 
dynamical Weiss mean-field and is determined in a self- 
consistent manner. Here we use the Nambu notation 
b^(r) = (6^ ) (r),&^ ) (r)*). Moreover the static bosonic 
mean-fields are defined in terms of the bosonic operator 
bj i/ as 

4?» = ^)o- (3) 

The index indicates that all averages are taken for the 
cavity system i.e. excluding the impurity site. Now each 
of the local actions can be treated as an impurity in the 
presence of a bath (representing the influence of the rest 
of the lattice) and therefore captured via an Anderson 
impurity Hamiltonian. There are several techniques to 
solve the impurity model. Here we apply Exact Diago- 
nalization (ED) [261 ] . 

In practice, we start with an initial set of local Weiss 
Green's functions and local bosonic superfluid order pa- 
rameters (frflir). After solving the action , we obtain 

a set of local self-energies S^(zw„) with u) n being Mat- 
subara frequency. Then we employ the Dyson equation 
in real-space representation in order to compute the in- 
teracting lattice Green's function 

G^r 1 = Go(^„)" 1 - S(twn). (4) 
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The site-dependence of the Green's functions is shown 
by boldface quantities that denote a matrix form with 
site-indexed elements. Here Go(iw„) _1 stands for the 
non-interacting Green's function 

Goiiuin)- 1 = (/x + iw„)l-t- V. (5) 

In this expression, the matrix elements are hopping 
amplitudes for a given lattice structure and the external 
potential is included via Vij = SijVi with Vi = V$rf and 
Ti being the distance from the trap center. Eventually 
the self-consistency loop is closed by specifying the Weiss 
Green's function via the local Dyson equation 

(<L(^n)) _1 = (G^^y 1 + EgttWn). (6) 

where the diagonal elements of the lattice Green's 
function yield the interacting local Green's function 

G^l(iu> n ) — (Gx,iy{iuJn))ii- This self-consistency loop is 
repeated until the desired accuracy for superfluid order 
parameters and Weiss Green function is obtained. 

Complementary to RBDMFT, in this work we employ 
an LDA approach combined with single-site BDMFT to 
explore the physics of the model ([TJ). Here we adjust the 
chemical potential on each lattice site according to LDA 
as ^i/(r) = yP v — V{r). The advantage of this approach 
is the larger system size accessible. We validate it by 
comparison with the more rigorous RBDMFT approach. 
In this paper, we apply RBDMFT for 2D lattices and 
LDA+BDMFT both for 2D and 3D lattices. 



III. MAGNETIC PHASES OF A 
HOMOGENEOUS BOSE-HUBBARD MODEL 
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FIG. 1: Upper panel: zero-temperature phase diagram for 
the two-component Bose-Hubbard model in a 3D cubic lattice, 
as a function of hopping parameters. The interaction values 
are Ub = Ud — 1.01 Ubd and the total filling is n = 1 with 
nb — nd = 0.5 (except in the phase separation regime). Lower 
panel: finite temperature phase diagram (T = O.OOlUbd). The 
energy scale of the magnetic coupling 4tbtd/Ubd is shown by 
the black dashed line. 



We start by exploring the two-component Bose- 
Hubbard model in the homogeneous case. We consider a 
3D cubic lattice and focus on the situation of total filling 
n = rib + n d being n = 1 and n = 2 with balanced den- 
sities rib — fid — 0.5 and rib = nd = 1 respectively. For 
each filling, we calculate both zero and finite tempera- 
ture phase diagrams. We focus on the interaction regime 
where the inter-species interactions, Ubb = Ub, Udd = Ud, 
and intra-species interaction, Ubd, are just slightly dif- 
ferent, i.e. Ub = U c i = 1.01 Ubd- This particular regime 
of interactions is accessible by Feshbach resonances, and 
indeed our choice is motivated by a recent experiment 
at MIT [13] , where a sample of 87 Rb atoms in hyperfine 
states |1, — 1) and |2, — 2) with nearly equal inter- and 
intra-species interactions has been prepared. The selec- 
tion of interactions {Ubd slightly smaller than Ub.d) is due 
to more novel magnetic phases appearing in this regime 
In all our calculations we set Ubd = 1 as the unit of 
energy, and z is the number of the nearest neighbors for 
each site. 

Fig. Q] displays the zero and finite temperature phase 
diagrams of the system with a total filling of one particle 
per site, n = 1. At T = (upper panel) we find three 



distinct phases which are characterized according to the 
value of the superfluid order parameters <f>b, 4>d arid the 
two-body correlator <pbd — (b<P) — (b)(cft) > which in- 
dicates the XY-ferromagnetic spin-ordering. When both 
species have comparably large hopping, we find a super- 
fluid phase characterized by (f>b,d 7^ 0. Instead when the 
hopping amplitudes are very different, the species with 
larger hopping is more easily delocalized and therefore 
superfluid, while the other component favors a Mott in- 
sulating phase. In this parameter regime, we do not 
find a homogeneous converged BDMFT solution where 
each component has the same filling, which indicates 
phase separation between the superfluid and the Mott 
insulator. We notice that the phase diagram is sym- 
metric upon particle interchange and this symmetry is 
also manifested in the Hamiltonian ([T]). The third phase 
emerges when the hopping amplitudes are small. This 
non-superfluid (i.e Mott insulating) phase possesses an 
XY-ferromagnetic spin-ordering and is characterized by 
4>b,d = and 4>bd > 0. 

We investigate also the effect of finite temperature on 
the phase diagram as shown in the lower panel of Fig. [TJ 
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FIG. 2: Upper panel: zero-temperature phase diagram for a 
two-component Bose-Hubbard model in a 3D cubic lattice, as 
a function of hopping parameters. The diagonal dotted line 
shows t b = td- The interaction values are U b = Ud = 1.01 U b d 
and the total filling is 2 with n b — rid = 1. Lower panel: finite 
temperature phase diagram (T = 0.0035U b d). 



We observe that the superflmd remains robust against 
small finite T. On the other hand, the XY-ferromagnetic 
spin-ordered phase is sensitive to finite temperature since 
it is formed in the low-hopping regime and therefore eas- 
ily destroyed by thermal fluctuations. At finite T this 
ordered phase is reduced in favor of developing a non- 
magnetically ordered Mott state which is characterized 
by <f>bd = and 4>b,d = and which we denote as "un- 
ordered insulator" in the following. The black dashed line 
shows the energy scale of the magnetic coupling ^htd/Ubd 

Now we turn to the case when the total filling at each 
site is n = 2. Fig. [5] (upper panel) shows the zero- 
temperature phase diagram for this case. The main dif- 
ference compared to n = 1 is the presence of a large 
unordered (Mott) insulator at low hopping values. As 
for n = 1, here we also find a sizable superfluid regime 
when both species have large hopping amplitudes. In- 
stead, when the hopping amplitude for one component is 
very small and the other one very large, the system will 
be in a new phase with one component being superfluid 
and the other one Mott insulating. This one- component 
superfluid phase, e.g. for the d component, is defined 
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FIG. 3: Zero-temperature phase diagram for a two- 
component Bose-Hubbard model in a 3D cubic lattice, as 
a function of hopping parameters. The interaction values 
are U b — Ud — 12 U b d and the total filling is n = 1 with 
n b = nd — 0.5. 



by 4>b = and 4>d > 0. There are also two Mott states: 
the XY-ferromagnet at intermediate hopping and the un- 
ordered (Mott) insulator in the lower hopping regime. At 
finite temperature (lower panel) both superfluid phases 
remain robust and almost unchanged. The main effect of 
finite T is to reduce the XY-ferromagnetic phase in favor 
of the unordered Mott insulator. The finite- T phase di- 
agram of Fig. [2] is in remarkable agreement with the one 
obtained by a field-theoretical approach [ilj . 

A further important spin-ordered state is the anti- 
ferromagnetic (AF) phase. Its existence in the Mott do- 
main has been shown in previous investigations (Til [T^j . 
In order to address this phase within our formalism, we 
adopt a set of parameters in which the inter-species inter- 
actions Ub,d are much larger than the intra-species one: 
Ub = Ud = 12 Ubd- This specific choice is inspired by 
a previous BDMFT study (26| in which the phase dia- 
gram was obtained on the Bethe lattice. Here we aim 
to map out the phase diagram on a 3D cubic lattice 
which is directly relevant for the experimental studies. 
Fig. [3] sketches the phase diagram for this case with the 
total particle filling n = 1. In addition to two previ- 
ously discussed phases, superfluid and XY-ferromagnet, 
(see Fig. (TJ [5]), we find two other ordered states here: 
AF phase and super-solid. For unequal hopping in the 
Mott domain, we identify a magnetically ordered phase 
of AF type. This non-superfluid phase (i.e. <pb,d = 0) is 
characterized by a finite value of the AF order parameter 
^af = \ n v,a — n v .a\ > 0, where v denotes the component 
and a is the sublattice (a = —a), together with vanishing 
XY-ferromagnetic order 4>bd = 0. Finally for a very large 
difference in the hopping of the two species we observe a 
small window of a super-solid phase with <j>b > 0, 4>d — 
and Aaf > if t b 3> td, and vice versa. 

To investigate in detail the quantum phase transition 
into the XY-ferromagnetically ordered state, for n = 2 
and at T — 0, in Fig. 0] we plot the behavior of individ- 
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FIG. 4: td dependence of order parameters (b), (d), and 
along the diagonal black dotted line in the upper panel of Fig. 
[2] The interaction regime is set to Ub = Ud = 1.01 Ubd and the 
hopping amplitudes it, = td with filling factors nt — rid = 1. 



ual superfluid order parameters and also the correlator 
</>6rf along the line of 4 = td (the diagonal black dotted 
line shown in the upper panel of Fig. [5]). The latter 
behavior indicates a second-order quantum phase transi- 
tion from the Mott state to the XY-ferromagnet and also 
a second-order phase transition from XY-ferromagnet to 
superfluid. 

One crucial question regarding observation of the AF 
and XY-ferromagnetic orders is how fragile they are 
against finite-temperature effects. To address this im- 
portant issue, we compute the respective critical tem- 
peratures. Fig. [5] shows T c as a function of the hopping 
amplitude of species b while we keep the hopping ratio 
fixed as t d = it b for the AF phase (upper panel) and 
as tb — td for the XY-ferromagnetic phase (lower panel) . 
We notice that T c rises as the hopping amplitudes in- 
crease, due to the growing effective exchange couplings 
which stabilize magnetic order. We also note that the 
zero-temperature phase diagram on a 3D cubic lattice for 
filling n = 1, containing the AF phase (Fig. [3]), has the 
same structure as the corresponding one on the Bethe 
lattice [26|. Therefore we anticipate that the finite- T 
counterpart of this phase diagram should also be similar 
on both lattices, and therefore expect there is a region 
of unordered Mott insulator at low hopping. The inset 
of Fig. [5] (upper panel) shows the temperature depen- 
dence of A^ F - It indicates a second-order phase transi- 
tion from the AF phase to the unordered Mott insulator. 
We have also computed the order parameter cftbd for the 
XY-ferromagnetic phase as shown in the lower inset of 
Fig. [5] which indicates a 2nd order transition from the 
XY-ferromagnetic phase to an unordered Mott insulator 
as well. The critical temperatures of magnetic phases 
shown here are notably smaller than the coldest temper- 
atures which have been measured in most experiments 
until now, apart from W. Ketterle's group where temper- 
atures as low as 350 pK (w O.OlUbd) have been achieved 
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FIG. 5: Critical temperature of AF and XY-ferromagnetic or- 
der as a function of the hopping amplitude tb on a 3D homo- 
geneous cubic lattice with total filling n — 1. Upper panel: 
AF phase with hopping amplitude ratio tb = 4t<j. Inset: melt- 
ing of the AF phase vs temperature with hopping amplitudes 
2ztb — O.QUbd and 2ztd = 0.15Ubd- Lower panel: XY- 
ferromagnetic phase for equal hopping amplitudes tb = td- 
Inset: melting of the XY-ferromagnetic phase vs temperature 
with hopping amplitudes 2ztt = 2ztd = 0.1875 Ubd- 



IV. SPIN-ORDERING FOR AN 
INHOMOGENEOUS BOSE-HUBBARD MODEL 

In the previous section, we focused on homogeneous 
systems. However all the experiments are carried out in 
the presence of an external confining potential. Therefore 
we extend the BDMFT scheme to real-space BDMFT to 
address the inhomogeneous system. In this section, using 
RBDMFT and LDA+BDMFT, we will explore an inho- 
mogeneous two-component Bose-Hubbard model both in 
2D and 3D. 



A. 2D trapped Bose gas 

In this section we discuss the AF phase, the XY- 
ferromagnet and the unordered Mott state in a 2D square 
lattice in the presence of a harmonic trap. We first in- 
vestigate AF ordering on a 31 x 31 lattice and then XY- 
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FIG. 6: Neel-type AF order in 2D at central filling n = 1 
and T = 0. Top: particle densities of the two species as 
a function of radial distance r. The interactions are set to 
Ub = Ud = 12(7^, the hopping amplitudes 2zt b — O.lUbd 
and 2ztd = 0.25 17m with harmonic trap Vb = 0.01 Ubd- The 
chemical potentials are /if, = /i^ = 0.5 ?7f,d- Bottom: density 
distribution of the d component over the lattice. 



FIG. 7: Neel-type AF order in 2D at central filling n = 2 
and T = 0. Top: particle densities of the two species versus 
radius obtained by RBDMFT. Interactions are set to Ub = 
Ud = 48 Ubd, the hopping amplitudes to 2ztb = 0.1 Ubd and 
2zf<i = 0.25 U b d, with a harmonic trap strength Vb = 0.01 Ubd- 
The chemical potentials are fib = — 1-5 ^m- Bottom: 
density distribution of the d component over the lattice. 



ferromagnetic phase and the unordered Mott state on a 
32 x 32 lattice. The choice of different lattice sizes is 
solely due to computational convenience. 



1. Anti-ferromagnet order 

One of the most desirable goals in the current exper- 
iments on cold atomic-gases is to reach the regime of 
(Neel-type) AF ordering which is (for fermionic systems) 
expected to be a key step towards realizing a d-w&ve su- 
perfluid phase @. Here we investigate this phase for 
trapped two component bosons in an optical lattice. At 
the beginning we focus on AF order at T = 0. We choose 
the interaction parameters as Ub = Ud = 12 Ubd with the 
unequal hopping for two species fixed as 2ztb = 0.1 Ubd 
and 2ztd = 0.25 Ubd- We choose a maximum local filling 
of n = 1 at the center of the trap. The top panel of 
Fig. [5] shows the RBDMFT results for the particle den- 
sities as a function of radial distance r from the trap 
center. The AF phase forms in the central area of the 
lattice as a checker-board pattern and vanishes smoothly 
with increasing distance r from the lattice center. This 
indicates that AF order is stable at the center of a finite 
trap. However due to the unequal hopping, the lighter 
species (i.e. the one with the larger hopping) explores the 
lattice more freely and forms a superfluid ring around the 



central checker-board pattern. This behavior is visible in 
the bottom panel of Fig. [5] To see how robust AF order 
is against changing the total atom filling number, we now 
increase the filling at the center of the trap to n = 2, and 
observe that AF order now forms as a ring around the 
center as shown in Fig. \7\ We conclude that also in 2D 
and in the presence of a trapping potential, AF order can 
be observed in regions of total filling n = 1 at T = 0. 



2. XY-ferromagnet 

As evident in the phase diagrams of the homogeneous 
system (Figs. [TJ [2]), a common magnetic phase which 
appears for both fillings n = 1 and n = 2 is the XY- 
ferromagnet. Here we investigate the stability of the XY- 
ferromagnetic phase in a trapped 2-component system on 
a 2D lattice of size 32 x 32 at T — 0. We first focus on the 
case of equal hopping 2ztb — 2ztd = 0.175Ubd for both 
species and choose the interactions as Ub — Ud — l.OlUbd- 
In Fig. [5] the atom densities, their corresponding su- 
perfluid order parameters and the correlator (f>bd are 
shown as determined from RBDMFT (top panel) and 
LDA+BDMFT (lower panel). At the center of the lat- 
tice, we have a total filling of n = 2. We observe a finite 
value of the correlator 4>bd which implies a stable XY- 
ferromagnetic phase in this inhomogeneous system. With 



7 




r 



FIG. 8: Top: RBDMFT results in 2D for atom densities, 
supcrfluid order parameters and XY-ferromagnetic correlator 
4>bd as a function of radial distance r at T = 0. Interac- 
tions are set to Ub = Ud = 1.01 Ubd, with hopping amplitudes 
22t(, = 2ztd = 0.175 (7td and harmonic trap Vo = 0.01 Ubd- 
The chemical potentials are fit — fid = 1.5 C/bd. Bottom: 
LDA+BDMFT results for the same parameters as the top 
panel. 
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FIG. 9: Unordered Mott state in 2D at n = 2 for weak hop- 
ping amplitudes at T = 0, calculated by RBDMFT (top 
panel) and LDA+BDMFT (bottom panel). The interac- 
tions and hopping amplitudes are Ub = Ud = 1.01 Ubd and 
Iztb — 2ztd — 0.1 Ubd, with the harmonic trap Vo = 0.01 Ubd- 
The chemical potentials are Ub — fid = 1.5 Ubd- 



increasing distance from the trap center, we find non-zero 
values for the superfluid order parameters with a maxi- 
mum inside the atomic cloud, indicating the superfluid- 
ity for both species. We can also see that the correlator 
<pbd remains finite in the superfluid regime. Moving fur- 
ther towards the edge of the trap, the XY-ferromagnetic 
phase with n = 1 appears and eventually a further su- 
perfluid domain. For comparison, we have computed the 
same quantities within LDA+BDMFT, which are shown 
in the bottom panel of Fig. [H We find a excellent agree- 
ment between RBDMFT and LDA+BDMFT deep inside 
each phase, but RBDMFT provides the more accurate 
description of the smooth transition between the differ- 
ent phases. Note that the second superfluid ring is very 
narrow within RBDMFT which is most likely a finite-size 
effect. 



3. Unordered Mott state 

We now consider the case of low hopping amplitudes. 
This situation corresponds to the regime of the large area 
of the unordered Mott state (without symmetry break- 
ing) in the homogeneous phase diagram (see Figs. lll2l for 
the 3D case) . In Fig. [5] we show results for the atomic 
densities, superfluid order parameters, and the correla- 



tor (pbd for the case 2ztb — 2ztd = 0.1 Ubd at T = 0. 
In the center of the lattice where we have filling n = 2, 
the low-energy Hilbert space of each lattice site includes 
the three states \b,b),\d,d) and \b, d). With the choice 
of the interactions as Ub,d > Ubd, the third state has the 
lowest energy and therefore we have a Mott state with 
4>bd = and no symmetry breaking. For intermediate 
radii, where 1 < n < 2 and both species are superfluid, 
(j>bd rises to a finite value and remains constant in the 
n = 1 region where the XY-ferromagnet is stable. 



B. 3D trapped case 

In this final part, we consider the experimentally most 
interesting case of a 3D cubic optical lattice in the pres- 
ence of an external harmonic trap. Due to the compu- 
tational limitations for RBDMFT, here we only apply 
LDA+BDMFT which we previously benchmarked ver- 
sus RBDMFT. Throughout this section, we consider a 
lattice with 41 x 41 x 41 sites. 
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FIG. 10: Top: particle densities, superfluid order param- 
eters and XY-ferromagnetic correlator as a function of ra- 
dius r on a 3D cubic lattice at T = 0, calculated within 
LDA+BDMFT. Interactions are set to U b = U d = 1.01 U bd , 
hopping amplitudes 2zt b = 2zt d = 0.195 Ubd with a harmonic 
trap V = 0.005 U bd and N tot = 2.6 x 10 4 . Bottom: Temper- 
ature dependence of the XY-ferromagnetic correlator for the 
same parameters as the top panel. 



1. XY-ferromagnet 

We begin the 3D trapped lattice analysis by investigat- 
ing the stability of the XY-ferromagnetic phase. We first 
choose intermediate hopping as 1zt\, = 2ztd = 0.195 Ubd- 
This choice corresponds to the XY-ferromagnet in the 
homogeneous phase diagrams (Fig. 1 1121) . We enforce fill- 
ing n = 2 at the center of the lattice by adopting the 
proper chemical potentials. In Fig. [10] (top panel) we 
show the particle densities, superfluid order parameters 
and the correlator fad at T — (top panel) and at finite 
T (bottom panel). For T = we observe a wedding- 
cake structure of the atomic densities i. e. two plateaus 
of n = 1 and n = 2 and intermediate areas with non- 
integer filling. However we are more interested in the 
magnetic properties which are revealed by the correlator 
fad- We observe that in the n = 2 domain, there is a 
XY-ferromagnetic phase, manifested by a finite value of 
fad and vanishing superfluid order parameters. In the in- 
termediate area 1 < n < 2, we obtain a superfluid phase 
with both (b) and (d) being finite. Note that the onset of 
superfluidity leads to non-zero XY-ferromagnetic corre- 
lations as well. By approaching the second Mott plateau 
with n = 1, the superfluid order parameters vanish again 



and we obtain a non-zero value for fad, indicating once 
again an insulating XY-ferromagnet. Finally for n < 1 
we find a further superfluid domain. 

As in the previous section, we are interested in the ef- 
fect of temperature on magnetic order. Fig. [TO] (bottom 
panel) represents the correlator fad for different temper- 
atures. First we notice that for all the temperatures 
considered here, the correlator possesses a larger value 
at n = 2 compared to n = 1. In other words, XY- 
ferromagnetic order is more pronounced for n = 2 com- 
pared to rt = 1 as long as all other parameters of the Bosc- 
Hubbard model are identical. To make this point more 
clear, we calculate the critical temperature for ferromag- 
netic order for both n = 1 and n = 2, and find respec- 
tively T c = 0.0018 Ubd (~ 70 pK) and T c = 0.0051 U bd 
(« 190 nK). We also calculate the maximum value of 
the critical temperature for XY-ferromagnetic order at 
filling n — 1 and n = 2, and find that the maximum 
value of the critical temperature is around 200 pK and 
230 pK, respectively, when the 3D cubic lattice is formed 
by laser beams of wave-length 1064 nm and the scatter- 
ing length is around 100ab (a b is the Bohr radius). This 
fact could be significant for ongoing experiments, e.g. in 
Refs. [22, 23[ where spin gradient thermometry has been 
used to measure temperatures as low as 350 pK in a 3D 
optical lattice. Our calculation here indicates that it is 
much easier to observe XY-ferromagnetism for higher fill- 
ing due to the enhanced critical temperature. 
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FIG. 11: Particle densities and XY-ferromagnetic correla- 
tor vs radial distance r for asymmetric hopping amplitudes 
2zt b = 0.15 U b d and 2zt d ~ 0.225 U bd on a 3D cubic lattice. 
Interactions are U b = U d — 1.01 U bd , the trapping potential is 
V = 0.005 Ubd and N toi = 2.6 x 10 4 . 

We now consider asymmetric values for the hopping 
amplitudes 2ztb = 0.15 Ubd and 2ztd = 0.225 Ubd- By ad- 
justing the chemical potentials, we obtain a globally al- 
most balanced mixture with N b ~ N d ~ (48%-52%) N tot 
with Ntot — 2.6 x 10 4 . Fig. [TT1 shows the atomic densities 
and XY-ferromagnetic spin-order for different tempera- 
tures (four panels). The most remarkable new feature 
of the asymmetric hopping regime is the vanishing ferro- 
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FIG. 12: Z-magnetization (imbalance) S z vs temperature for 
the homogeneous system. Interactions are set to Ut = Ud — 
1.01 Ubd with hopping amplitudes 2ztt = 0.15 Ubd and 2ztd = 
0.225 U bd and total filling n = 2. 



magnetic order in the n = 1 domain where = even 
at T = 0. We also find that the total density profile first 
becomes sharper and then smoother again with increase 
of temperature. This is due to the higher spin entropy 
of the unordered two-component Mott insulator. Details 
will be discussed in a future publication [34| . 

It is visible in Fig.[TT]that the lighter species (i.e. the 
one with larger hopping) always dominates the density 
distribution in the superfluid area where n ^ 1,2, since in 
this regime the particle mobility plays a more important 
role than in the Mott domains. For the Mott state (n = 
1, 2) the situation differs. In the spin model language 
\12\ , different hopping amplitudes lead to the existence of 
an effective magnetic field [13, E3 , 

h = z(2S+l)^— ^+ fi b ~ fi d . 

Ubd 

where S = n/2. This effective magnetic field gives rise 
to an imbalance between rib and rid which we quantify 
by the Z-magnetization S z = (rib — "d)/2. At the Mott 
plateau with n = 1, the magnetization will shrink from 
1/2 (maximum of Z-magnetization due to a large effective 
magnetic field) to zero with increasing T due to thermal 
fluctuations. At the Mott plateau with n = 2, the mag- 
netization depends non-monotonically on temperature. 
In Fig. [T21 we show the temperature dependence of the 
magnetization at filling n = 2 plateau by focusing on 
the trap center and performing a finite-T study with a 
single-site BDMFT for the homogeneous model. This 
behavior can be understood if we notice that in the fill- 
ing n = 2 region the system favors a state with a small 
magnetization S z at zero temperature due to a non-zero 
effective magnetic field. Once the temperature starts to 
increase from zero to a finite value, thermal fluctuations 
will come into play and compete with quantum fluctua- 
tions, which makes the imbalance reach a maximum at 
finite T, since both types of fluctuations can delocalize 
the atoms and smooth the imbalance between the two 
species. When the temperature increases even higher, 



the larger thermal fluctuations will simply smear out the 
imbalance. 



2. Unordered Mott state 




r 



FIG. 13: Unordered Mott state at n = 2 for weak hop- 
ping amplitudes at T = 0, calculated within LDA+BDMFT. 
Interactions are Ut = Ud = 1.01 Ubd, hopping amplitudes 
2ztb — 2ztd = 0.15 Ubd, with a harmonic trap Vb = 0.005 Ubd 
and Ntot = 2.6 x 10 4 . 



Finally, we also investigate the non-magnetic Mott 
state without symmetry breaking which occurs at n = 2 
for relatively low hopping amplitudes 2ztb — 2ztd — 
0.15Ubd (see also the homogeneous phase diagrams 
Fig. [2]). Fig. [13] (top panel) shows the results for atomic 
densities, their superfluid order parameters and the cor- 
relator <pbd- At the center of the trap we indeed find a 
Mott state without symmetry breaking similar to the 2D 
case (Fig. EJ. 



V. SUMMARY AND OUTLOOK 

In conclusion, we have studied magnetic ordering of a 
two-component Bose gas in 2D and 3D optical lattices. 
By using BDMFT we have calculated the phase diagrams 
of the homogeneous Bose- Hubbard model in a 3D cubic 
lattice with total particle filling n = 1 and n = 2, which 
feature superfluid and Mott-insulating phases and also 
notably reveal ordered phases with XY-ferromagnetism 
and Neel-antiferromagnetism in the Mott domain. We 
investigate the critical temperatures of these long-range 
ordered states. Moreover we have confirmed the sta- 
bility of these magnetic phases in a trapped 2D or 3D 
system. In the case of a 3D cubic lattice, we have in 
particular computed XY-ferromagnetic ordering at finite 
temperatures which is relevant for current experiments 
(22l [23J]. Another important issue is the detection of 
novel magnetic phases with long-range spin order. For an 
antiferromagnetic phase, spin-sensitive detection can be 
used to detect the Neel-type order, i.e. one spin compo- 
nent can be detected after removing the other one using 
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spin-selective single-site addressing, as pointed recently 
[111, HU- For the XY- ferromagnetic phase, one could ap- 
ply a 7r/2 pulse with a position-dependent phase (lin- 
ear gradient), thus probing locally with different phases. 
For XY-ferromagnetic long-range order one would thus 
observe stripes with regular spacing when the phase 
matches and the tt/2 pulse transfers all atoms into one 
of the two spin states [37|. 

Achieving the necessary ultra-low temperatures for de- 
tecting magnetic ordering of cold bosons in optical lat- 
tices has so far remained elusive. However, it is antici- 
pated that by further experimental advances this obsta- 
cle will be overcome in the near future 1381. Our results 



provide theoretical benchmarks which in the future will 
also be extended to other geometries such as triangular 
or hexagonal lattices. 
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